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We present the result of a search for the charmless two-body baryonic decay B" ^ pp in a sample 
of 88 million T(4S) BB decays collected by the BaBar detector at the SLAC PEP-II asymmetric- 
energy B Factory. We use Cherenkov radiation to identify protons cleanly, and determine the signal 
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yield with a maximum-likelihood fit technique using kinematic and topological information. We find 
no evidence for a signal and place a 90% confidence-level upper limit of B{B° ~* pp) < 2.7 x 10"^ 



PACS numbers: 13.25.Hw, ll.30.Er, 12.15.Hh 

We report the result of a search for the charmless two- 
body baryonic decay ^ pp Q- Although B mesons 
have recently been observed to decay into several charm- 
less three-body baryonic final states jH , there is currently 
no evidence for the corresponding charmless two-body 
decays. Previous searches 0, Q for B° pp decays 
have yielded upper limits on the branching fraction at 
the level of 10^^, which is consistent with calculations 
based on QCD sum rules and the pole model [^. 
A simple scaling of the measured branching fraction for 
B" A-p 13 by the current estimate [3 of |Kb/Vcb|^ 
leads to a prediction of charmless two-body branching 
fractions at the level of 10~^, which is near the current 
sensitivity of present experiments. 

The data sample used for this search contains (87.9 ± 
1.0) X 10^ T{AS) -> BB decays collected by the BABAR 
detector 9j at the SLAC PEP-II e+e^ asymmetric- 
energy storage ring. The primary detector components 
used in the analysis are a charged-particle tracking sys- 
tem consisting of a five-layer silicon vertex detector and a 
40-layer drift chamber surrounded by a 1.5-T solenoidal 
magnet, and a dedicated particle identification system 
consisting of a detector of internally reflected Cherenkov 
light (DIRC). 

Two-body B decays are reconstructed from pairs 
of oppositely-charged tracks originating from the in- 
teraction region and having momentum greater than 
lOOMeV/c in the direction transverse to the beam line. 
We require each track to have an associated Cherenkov 
angle (Oc) measurement with at least four signal photons 
detected in the DIRC. To suppress combinatorial back- 
ground arising from A decays, we require that the two 
tracks form a vertex with probability greater than 10"^. 

Signal candidates are identified kinematically with two 
variables: the difference AE between the center-of-mass 
(CM) energy of the B candidate and \/s/2, where y/s is 
the total CM energy, and the beam-energy substituted 
mass rriEs = \/ (s/S + Pi • y>bY /Ef — p^, with the B- 
candidate momentum pb and the four-momentum of the 
initial state {Ei, pi) defined in the laboratory frame. For 
signal decays, AE peaks near zero with a resolution of 
about 23 MeV, while toes peaks near the B mass with a 
resolutionof about 2.6 MeV/c^. We require 5.20 < toes < 
5.29GeV/c2 and |A£;| < 100 MeV. 

Protons are identified based on the 6^ measurement 
from the DIRC and the momentum measurement from 
the tracking system. Figure shows the difference be- 
tween measured and expected values of 9c for the pro- 
ton hypothesis, divided by the error ug^, for tracks from 
B'^ ^ pp candidates in the sideband region 5.20 < toes < 
5.26GeV/c^. Protons peak near zero and are well sep- 
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[e^-e^(exp)]/ag^ 

FIG. 1: The difference between the measured and expected 
values of Oc, divided by the error, for tracks from _B° — > pp 
candidates in the region 5.20 < mEs < 5.26GeV/c^ We only 
use tracks that lie on the left side of the dashed line. 



arated from the much larger background of pions and 
kaons. We require a 9c measurement within hag^ of the 
expected value for a proton, which removes over 97% of 
the combinatorial background while retaining more than 
91% of the signal decays (the efficiency is less than 100% 
due to the presence of non-gaussian tails in the pull dis- 
tribution). 

We measure the efficiency of the 9c selection in a sam- 
ple oi A^ pTT~ decays reconstructed in 9.6 fb~^ of e+e" 
annihilation data recorded 40 MeV below the T{AS) reso- 
nance. The sample is selected using kinematic and decay- 
vertex information, and has a purity of 98.5%. For con- 
sistency with B^ pp decays, we require the proton CM 
momentum p* to be in the range 2.2 < p* < 2.8 GeV/c. 

Due to the unique topology and kinematics of the two- 
body final state, b c decays do not populate the signal 
region for —^ pp, and backgrounds from b ^ u de- 
cays are negligible after the proton selection. We verify 
both assertions by analyzing a sample of approximately 
80 fb"^ of 7(45") BB Monte Carlo simulated events in 
which the B mesons decay according to the world-average 
branching fractions 01 , and a second sample correspond- 
ing to approximately 200 fb~^ where one B meson in each 
event is forced to decay to a charmless final state. No 
event passes the above selection requirements in either 
sample. 

The dominant background is from random combina- 
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tions of protons produced in the process e+e~ qq 
{q = u,d,s). We verify in Monte Carlo samples that 
the background from the process e+e~ cc is negligible 
compared to light-quark production. In contrast to the 
spherical topology of BB events, particles produced in 
light-quark events tend to lie near the thrust axis of the 
original qq pair. To suppress this background, we cal- 
culate the angle 9s between the sphericity axis of the B 
candidate and the sphericity axis of the remaining par- 
ticles in the event, and require | cos 6*5 1 < 0.9. This re- 
quirement removes 70% of the combinatorial background, 
while retaining 85% of the signal decays. In addition, we 
define a Fisher discriminant |lOl |. which is a sum of 
two discriminating variables with coefficients optimized 
to separate signal and light-quark events. The first vari- 
able is the scalar sum of the CM momenta of all the 
particles in the event, excluding the two tracks from the 
B'^ pp candidate. The second variable is the prod- 
uct p* (cos 6**)^ summed over all particles (excluding the 
B-candidate tracks), where 9* is the angle between its 
momentum and the -B-candidate thrust axis in the CM 
frame. 

The total efficiency for all of the above selection criteria 
is (35.8 ± 3.7)%, where the error includes the statistical 
and systematic uncertainties added in quadrature. The 
dominant source of the uncertainty is due to the limited 
statistical precision of the A control sample after applying 
the proton p* constraint. A total of 804 events satisfy the 
B'^ — > pp selection criteria. 

The signal yield is determined from a maximum like- 
lihood fit that uses mss, and as discriminating 
variables. The likelihood for the sample is defined as 



N 



c 



(1) 



i=l 



where N is the total number of events in the sample, A''s 
and Nq are the signal (S) and background (B) yields, 
and and are the signal and background probabil- 
ity density functions (PDFs) evaluated for event i. The 
PDFs are calculated from the product of PDFs for the 
individual variables, which are taken to be uncorrelated 
in the fit. We verify this assumption by calculating the 
linear correlation coefficient for each pair of variables. 
The largest correlation (—13%) is between toes and AE 
in signal decays, and we have confirmed that the effect 
of this small correlation is negligible. The signal yield 
is determined by minimizing the function — 21n£ with 
respect to Ns and Nb- 

We use data and Monte Carlo samples to model the 
PDF shapes for signal decays. The mean and resolu- 
tion of Toes are dominated by the beam energy, and 
are therefore similar in decay modes where the momen- 
tum resolution of the B candidate is significantly better 
than the resolution on the beam energy. We obtain the 
mean and resolution of TOeS: and also the mean of AE, 



from a large sample of B~ D^{K~'k'^)-k~ decays re- 
constructed in data. Due to the difference in momen- 
tum resolution between protons and mesons, the reso- 
lution on AE is different for pp and D'^tt~ decays. We 
therefore use the value obtained in a large Monte Carlo 
sample of B° pp decays, and apply a 5% correction 
to account for the observed difference in AE resolution 
for D'^7t~ decays reconstructed in data and Monte Carlo 
samples. For we use an asymmetric gaussian function 
with parameters obtained from simulated events. The 
shapes of the background PDFs are obtained from data 
in the sideband regions 100 < |A£:| < 200 MeV and 
5.20 < Toes < 5.26GeV/c^. We use a linear shape for 
AE, a double-gaussian function for J^, and an empirical 
threshold function for toes 0|- 

Several cross-checks are performed to validate the fit- 
ting technique. To confirm that the signal yield is 
unbiased, we generate and fit a large set of pseudo- 
experiments where signal and background events are gen- 
erated randomly from the PDFs. For these studies, we 
assume a branching fraction of lO^*' and find that the 
fitted signal yield is unbiased. We also check for bi- 
ases arising from kinematic correlations by mixing Monte 
Carlo signal events with backgrounds generated directly 
from the PDFs. No significant biases are observed. The 
sensitivity of the analysis is determined from a set of 
pseudo-experiments with assumed branching fractions in 
the range (0.1-1.2) x 10~^. We find that for any branching 
fraction above 0.5 x 10~^, the null hypothesis would be 
excluded with a probability greater than 99.997%, corre- 
sponding to a significance of 5(t for a gaussian distribu- 
tion. 

The result of the fit is A's = — O.Slj'J' consistent with 
no signal. We determine a Bayesian 90% confidence-level 
(C.L.) upper limit on Ns by finding the value Ng^ such 
that 



Jo 



.dNs 



0.90, 



(2) 



where /!max is the value of the likelihood as a function of 
iVs- We find N^^ = 6.3 events. 

Figures[2Ia-c) show projections of the fit result in each 
of the discriminating variables. The data in the signal 
region agree well with the background PDF shapes de- 
termined from sideband data. As a cross-check on the 
fit, we apply more stringent background-rejection criteria 
and determine the signal yield from the observed num- 
ber of events in a restricted signal region in toes and AE'. 
We require |cos0s| < 0.7, and define the signal region as 
\AE\ < 30 MeV and 5.27 < toes < 5.29GeV/c2. Fig- 
ure HJd) shows the toes distribution for events passing 
the more restrictive cos 9s and AE requirements. There 
are 9 events in the signal region with an expected back- 
ground of 7.7 ± 1.4, where the background is determined 
by extrapolating the observed yield in the sideband re- 
gion 5.2 < Toes < 5.26GeV/c^. The signal yield is 
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FIG. 2: Data (points with errors) and the result of the maxi- 
mum likelihood fit (solid line) projected onto the (a) tties, (b) 
AE, and (c) J- variables, and (d) the distribution of tties for 
events in the signal-enhanced sample defined by the require- 
ments |cos0s| < 0.7 and 1A£;1 < 30MeV. 



TABLE L Summary of absolute systematic uncertainties on 
A'^s from variations in the PDF parameters. The total uncer- 
tainty is the sum in quadrature of the individual contribu- 
tions. 



Source Positive variation Negative variation 



Signal 






mES 


-1-0.11 


-0.62 


AE 


-fO.39 


-0.37 




4-0.03 


-0.03 


Background 






mES 


4-0.32 


-0.30 


AE 


-hO.Ol 


-0.01 


T 


-H0.86 


-0.85 


Total 


-H.OO 


-1.16 



TABLE IL Summary of relative statistical and systematic 
uncertainties on the signal efficiency. The total uncertainty is 
the sum in quadrature of the individual contributions. 



Source 


Uncertainty (%) 


Statistical 


7.7 


Tracking 


1.6 


Vertex Quality 


2.5 


Proton Selection 


3.0 


DIRC Acceptance 


1.0 


cos Qs 


5.0 



Total 10.2 



1.3 ± 3.3, which is consistent with the null result from 
the likelihood fit. 

Tables and ^ summarize the various sources of sys- 
tematic error on the signal yield and efficiency. System- 
atic uncertainty on No, may arise from imperfect knowl- 
edge of the PDF parameters. We vary each parameter 
by its estimated error and combine in quadrature the re- 
sulting variations in A^g . For the efficiency of the proton 
selection, we assign the 1.5% background fraction in the 
A control sample as the correlated systematic error. The 
efficiency of the vertex quality requirement is determined 
to be 97.5% from simulated pp decays, and we assign a 
systematic uncertainty of 2.5% to account for possible 
differences between data and Monte Carlo events. As a 
cross-check, we compare the efficiency in the topologi- 
cally similar decays B'^ ~* vr+vr^ and i?° K^tt^ and 
find good agreement between data and simulation. Fi- 
nally, we include a correlated systematic error of 0.8% 
per track to account for possible differences in tracking 
efficiency between data and Monte Carlo events. The to- 
tal systematic uncertainty on the efficiency is computed 
by adding correlated errors linearly, and then adding the 
separate sources in quadrature. 

We calculate the 90% C.L. upper Hmit on the branch- 
ing fraction by increasing iVg'^ by the total systematic 
error on the signal yield, and by decreasing the efficiency 
and number of BB events by their respective total uncer- 



tainties. We find the flavor- averaged branching fraction 
B{B° pp) < 2.7 X 10-^ at the 90% C.L. This result 
improves the previous limit by more than a factor of 
four. 

In summary, we have performed a search for the de- 
cay — > in a sample of 88 million BB events. We 
find no evidence for a signal and set an upper limit on the 
branching fraction at 2.7 x 10^^. This result rules out the 
calculation in |^ based on QCD sum rules, while it is con- 
sistent with a recent calculation using the pole model , 
and with simple scaling of the measured branching frac- 
tion for the decay B'^ — ^ ^cP- 
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